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Abstract 

We consider two parallel cyclic Ising chains counter-rotating at a rela- 
tive velocity v, the motion actually being a succession of discrete steps. 
There is an in-chain interaction between nearest-neighbor spins and 
a cross-chain interaction between instantaneously opposite spins. For 
velocities v > the system, subject to a suitable markovian dynam- 
ics at a temperature T, can reach only a nonequilibrium steady state 
(NESS). This system was introduced by Hucht et al., who showed that 
for = oo it undergoes a para- to ferromagnetic transition, essentially 
due to the fact that each chain exerts an effective field on the other one. 
The present study of the v = oo case determines the consequences of 
the fluctuations of this effective field when the system size N is finite. 
We show that whereas to leading order the system obeys detailed bal- 
ancing with respect to an effective time-independent Hamiltonian, the 
higher order finite-size corrections violate detailed balancing. Expres- 
sions are given to various orders in A^~^ for the interaction free energy 
between the chains, the spontaneous magnetization, the in-chain and 
cross-chain spin-spin correlations, and the spontaneous magnetization. 
It is shown how finite-size scaling functions may be derived explicitly. 
This study was motivated by recent work on a two-lane traffic problem 
in which a similar phase transition was found. 

Keywords: kinetic Ising model, nonequilibrium stationary state, 
phase transition 
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1 Introduction 



Recently Hucht |T] (see also [2]), motivated by the phenomenon of magnetic 
friction, formulated a nonequilibrium steady state (NESS) Ising model of a 
new type. It consists of two parallel linear Ising chains having a relative ve- 
locity V. In addition to a nearest-neighbor interaction in each chain, any pair 
of spins facing each other on the two chains has an instantaneous interaction. 
In the version of the model easiest to study, each chain is finite and periodic; 
we will therefore speak of cyclic counter-rotating Ising chains (CRIC). The 
model, subject to suitable temperature dependent Markovian dynamics, was 
shown at velocity f = cxd to have a para- to ferromagnetic phase transition 
which in the limit of infinitely long chains may be understood in terms of an 
equivalent equilibrium model. 

The CRIC seems to us to be of the same fundamental importance as 
Glauber's ^ original kinetic Ising model. First, it is of interest in its own 
right as a new member of the class of NESS. Second, its interest is enhanced 
in the wider context of recent work on Ising models that in one way or 
another are driven, dissipate energy, or have some novel type of coupling; 
such work has appeared in a variety of contexts [H El E] • In particular, the 
present CRIC was extended to a Potts version by Igloi et al. [7], who find 
remarkable nonequilibrium phase transitions. In this paper we contribute 
further to the study of the CRIC. We focus on finite chains and on how to 
derive known and new properties from the master equation that defines the 
model. 

Hucht 's solution [1] is based on showing that at f = oo the stationary 
state dynamics of the CRIC is actually that of an equilibrium Ising chain 
in an effective magnetic field Hq, this field being zero above the transition 
temperature and nonzero below. This equivalence is valid in the limit where 
the chain length tends to infinity. In this work we show that it is possible 
to formulate this problem as an expansion in powers of N"^/"^. To lowest 
order we recover the equivalent equilibrium system found in reference [1]. 
To higher orders fluctuations of the field Hq come into play and appear as 
finite-size effects. 

The finite case is of interest, first of all, on the level of principles, and 
secondly, for the analysis of finite size effects in simulations as were carried 
out in [1] and by ourselves. We expect, furthermore, that our approach will 
help prepare the way for future work on the v < oo case, which is considerably 
harder. 

The effective transition rates satisfy detailed balancing to leading order 
in the large- iV expansion [1]; our analysis reveals, however, that to higher 
orders in A^~^/^ the detailed balancing (DB) symmetry of the effective rates 
is broken. The stationary state distribution may be found explicitly, at least 
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Figure 1: Ladder of spins witli an intrachain nearest-neighbor interaction Ji = J. 
The two chains constituting the ladder have a relative velocity v, the motion taking 
place in discrete steps of one lattice unit. There is an interchain nearest-neighbor 
interaction J2 = r]J between each pair of spins facing each other at any instant in 
opposite chains. 

to the lowest DB-violating order. Knowing this state one can calculate all 
desired NESS properties. 

In section [2] of this paper we define the rules of the markovian dynamics 
for general relative velocity v and then specialize to f = 00. These dynamical 
equations are the starting point for all that follows. In section[3]we discuss the 
DB violation that occurs in higher orders of . In section H] we consider the 
stationary state to zeroth order, as was already done by Hucht p. In sections 
[5] and El we show how A^~^ can be introduced as an expansion parameter and 
we define a 'leading order', composed of the zeroth order and a first-order 
correction. In section [7] we show how for the stationary state distribution an 
expansion may be found in powers of A^~^. We present the explicit result 
to next-to-leading order. In section [S] we calculate for various quantities of 
physical interest their stationary state averages to successive orders in the 
expansion. In section [H] we briefly discuss the relation of the present model 
to a two-lane road traffic model studied earlier. In section [10] we conclude. 
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2 Counter-rotating Ising chains 



2.1 A stochastic dynamical system 

We consider Ising spins on the ladder lattice shown in figure [1] The spins 
in the upper chain are denoted by rj, those in the lower chain by Si, where 
the integers j and i are site indices. There is a nearest-neighbor interaction 
Ji = J inside each chain and an interaction J2 = rjj between each pair of 
spins facing each other in opposite chains. We take J > and rj of arbitrary 
sign. The feature [2] and [1] that distinguishes this model from the standard 
Ising model on a ladder lattice, is that the two chains move with respect to 
one another at a speed f > 0. This will mean the following: the time axis is 
discretized in intervals of duration r = ao/v (where oq is the lattice spacing) 
and at the end of each interval the upper chain is shifted one lattice spacing 
Oo to the right with respect to the lower one. The Hamiltonian 'H(t) of this 
system is therefore time-dependent and given by 

^(^) = -'^Yl [^J'^J+i + ^j^i+^l - ^i^l^+^tM ) (2-1) 

i j 

where \_x\ is the largest integer less than or equal to x. 

We will consider cyclic boundary condition^. In this case the chains 
become counter-rotating loops of length say A^; the site indices i, j, and 
[j + vt/ao\ must then be interpreted modulo A^. Employing the shorthand 
notation r = {rj\j = 1, 2, ... , A^} and s = {sj\j = 1,2,..., A^}, we may 
indicate a spin configuration of the system by (r, s). 

We associate with Hit) a stochastic time evolution of (r, s). Its precise 
definition requires that we exercise some caution. We will first define it as a 
Monte Carlo procedure and then write down the master equation and pass 
to analytic considerations. Single-spin reversals are attempted at uniformly 
distributed random instants of time at a rate of I/tq per sitd. Each attempt 
is governed by transition probabilities. Since there are 2A^ sites, there are 
2A^ different single-spin fiips by which a state (r, s) may be entered or ex- 
ited. Given that a reversal attempt takes place, let {2N)~^WJ{r] s;t) and 
{2N)~^Wj{s] r; t) be the probabilities that rj and Sj are flipped, respectively. 
The reversal attempt will remain unsuccessful with the complementary prob- 
ability 

N 

1 - A,ec = 1 - {2N)-' J2 [W-ir- s- 1) + W^{s- r- 1)] , (2.2) 

"'^In connection with the traffic problem open boundary conditions are certainly also 
worthy of consideration. These have however the inconvenience of breaking the transla- 
tional symmetry. 

^We may scale time such that tq = 1. 
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where A^^c is what is usually called the 'acceptance probability'. 

We now specify the Wj and Wj in such a way that at any time t the 
system strives to attain the canonical equilibrium at a given temperature 
T with respect to the instantaneous Hamiltonian Hit). The choice is not 
unique. We choose 

W]{r;s;t) = i [l - irj(rj_i + r^+i) tanh 2ir] [l - r^s^ tanh r/iT] , 

Wl{s;r;t) = i [l - isi(si_i + s^+i) tanh 27^] [l - s^r,- tanh r/iT] , (2.3) 

where we have set K = J/T (with T measured in units of Boltzmann's 
constant) and where in both equations i and j are related by 

i= [j + vt/ao\ mod A^. (2.4) 

Equation fl2.3p is different both from the heat bath (or: Glauber) and from 
the Metropolis transition probabilities. We will refer to it as the "factorizing 
rate" . The factor 

wf{r) = I [1 - ir,(rj_i + r,+i) tanh 2K] (2.5) 

represents the Glauber transition probability. The Wj and Wj define an easy- 
to-simulate Markov chaiij^l with time-dependent transition probabilities^ 

In the special case v = the Hamiltonian 'H(t) reduces to the equilibrium 
Hamiltonian of the ladder lattice. For v arbitrary but rj = it reduces to the 
equilibrium Hamiltonian of two decoupled chains. In both special cases the 
dynamics is standard and obeys detailed balancing. 

In the general case, since the Hamiltonian is time-dependent, the system 
will not reach equilibrium but instead enter a NESS. Actually, for generic 
V, because of the periodic discrete shifts, the NESS is a r-periodic function 
of time; NESS averages are naturally defined to include an average over this 
period. In the limiting case v = oo we have r = and this complication 
disappears. The infinite velocity NESS is the subject of our interest in the 
remaining sections. It is a problem that depends only on the two parameters 
K and rj. 

We note finally that as compared to ours, there is an extra prefactor 
2 + (l-r,_ir,+i)tanh2ir 

l + tanh2i^ ^^^"^ ' ^^-^^ 

in Hucht's expression for the transition probability Wj{r; s;t), and an anal- 
ogous prefactor for VTf (s; r; t). These factors may easily be carried along in 
the calculation. 



■'No confusion should arise with the two legs of the ladder lattice, to which we refer 
also as 'chains'. 

^The reversal attempts, that is, the steps of the Markov chain, are Poisson distributed 
on the time axis. This makes it possible at any time to probabilistically connect the 
elapsed time t to the number of spin reversal attempts n. In the large t limit of course 
n ~ t/TQ. 
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2.2 The limit v ^ oo 



Let P(r, s; n) be the probability distribution on the configurations (r, s) after 
n spin reversal attempts. We will now write down the formal evolution 
equation for P(r, s; n) for the case oiv = oo, where important simplifications 
occur. When v = oo there is no relation between the indices i and j and 
hence the chain has transition probabilities Wj{r] s) given by the average of 
(12.31) on all z, which is now considered as an independent variable. We denote 
this average by Wj(r, s) and thus have 



1 ^ 



i=l 



^ [l — \rj{rj_i + Tj+i) tanh2fC] [l — rjix{s) taiohr]K~\ 
wf{r) X i [l - rjfi{s) tanh'qK] , (2.7a) 



1 

N 



where 



N 

wf{s) X i [1 - Sj/i(r) tanhT^ir] , (2.7b) 
1 ^ 1 ^ 

1=1 i=l 

We will write for the configuration obtained from r by reversing rj (that is, 
by carrying out the replacement rj i— )■ —Vj), and similarly define . Summing 
on all 2N fiips by which it is possible to enter or to exit (r, s) we find that 
the evolution of P(r, s; n) is described by the master equation 



P(r, s;ri+l) = [ w.,(r-^; s)P{r^ , s; n) + Wj{s^] r)P{r, s^; n) 

i=i 

+ (l - Wj{r; s))P{r, s; n) + [l - Wj{s; r))P{r, s;n)] , 

(2.9) 

where the second line corresponds to the probability of an unsuccessful spin 
reversal attempt. In vector notation equation (12.91) may be written 

P(n+ 1) = (1 + W)P(n), (2.10) 

where P(n) is the 2^^-dimensional vector of elements P{r,s;n), the symbol 
1 denotes the unit matrix, and W is a matrix composed of entries wj for 
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which comparison of fl2.9p and fl2.10p yields 



W{r, S;r' , s') = 6j,irjSs'sU)j{r^ , s) + 5r'rSs's:>Wj{r, s^) 

N 

- Sr'rSs's ^ s) + Wj{s] r)] . (2.11) 

j=i 

The discrete-time master equation, (12. 9p together with the Poisson statistics 
of the reversal attempts on the time axis, fully defines the CRIC for v = oo. 
This equation may be studied analytically, as is the purpose of this work, or 
may be implemented in a Monte Carlo simulation. 

3 Detailed balancing and its violation 

Henceforth we consider the case f = oo. Our purpose is now to find the 
stationary state distribution Pst{r,s) of the evolution equation (12. 9p . This 
distribution is the solution of P{r,s;n) = P{r,s;n+ 1) = Pst{r,s), which 
means 

= WPst. (3.1) 

Combining equations (13. ip and (12.90 yields the v = oo stationary state equa- 
tion 

N 

= ^ [wj{r^; s)P^t{r\ s) + Wj{s^ ]r)P^t{r, s^) 

- Wj{r; s)Pstir, s) - Wj{s; r)Pst(r, s)] . (3.2) 

If the transition probabilities satisfy the condition of detailed balancing, the 
solution of (13.21) is easily constructed; in case of the contrary, there are no 
general methods. We examine therefore first the question of whether equation 
(12.91) satisfies detailed balancing. 

A Markov chain satisfies detailed balancing (DB) if and only if its tran- 
sition probabilities are such that any loop in configuration space is traversed 
with equal probability in either direction. To show that the transition proba- 
bilities Wj fail to obey DB we consider an elementary loop of four single-spin 
flips, 

(r, s) (r-', s) H-> (r-', s-') h-> (r, s^) t-)- (r, s). (3.3) 

Given the system is in (r, s), we denote by p+{ri) and P-{ri) the probability 
that in the next four attempts it goes through this loop in forward and in 
backward direction, respectively. That is, 

P+iv) = Wj{r; s)wj{s; r^)wj{r^; s^)wj{s^; r), 

P-iv) = Wj{s;r)wj{r;s^)wj{s^;r^)wj{r^;s). (3.4) 
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For 7] = the two chains are decoupled, and as discussed below equation 
(I2.4p . each of them separately satisfies DB; it is easy indeed to verify explicitly 
that p+{0) = P-{0) = p(0). For general rj we may work out the difference 
P+iv) ^ P-iv) with the aid of fl2.7al) . fl2.8p . and the relations 

2r • 2s- 

Mr^) = Mr)-^, Ks^) = Ks)-^, (3.5) 

which yields 

P+iv) " P-iv) = 4:N^^p{0) tanh^r^K [rjfi{s) — Sjfi{r)] 

x{[rjfi{r) - Sjfi{s)] + 2N-HanhriK} . (3.6) 

This shows that DB is violated in the general case of nonzero coupling (r/ 7^ 0) 
between the chains. It becomes valid again only asymptotically in the limit 
N ^ oc. We therefore cannot hope to rely on any general methods to con- 
struct Pst('^, s) for finite A^. Indeed, writing out the stationary state equation 
(13.11) fully explicitly for = 3,4 (only = 2 is trivial) has confirmed the 
nontriviality of the stationary state but has not provided us with any useful 
insight. 



4 Stationary state Pst{r,s) to zeroth order 

The limit N ^ 00 was considered by Hucht lU [2] , and we briefiy recall the 
results. One may suppose that in this limit n{r) and n{s) have vanishing fiuc- 
tuations around an as yet unknown common average to be called mo{K,ri). 
We will denote the N ^ 00 limit of wj by wj^ . It then follows from fl2.7al) 
that 

'^jfli^) = wfir) X |[1 — r^-mo tanhr^i^]. (4.1) 

With the transition probabilities (14. ip the r- and the s-chain decouple. More- 
over, the expression for these wj^ is such that the spin dynamics satisfies DB 
with respect to the pair of uncoupled nearest-neighbor Ising Hamiltonians in 
a field, 

AT AT 

Hoir, s)/T = -KJ2 h'^i+i + ^^•^.•+1] -HoJ2 h- + ^j] ■ (4-2) 

where Hq is defined in terms of mg by 

tanhifo = wiotanhr/A' (4.3) 

and where K and Hq both include a factor 1/T. The quantity TLo{r, s) is an 
effective time-independent Hamiltonian. Let m{K, z) denote the magnetiza- 
tion per spin of the one- dimensional (ID) Ising chain with coupling in a 
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field that we will for convenience denote by z. This quantity is well-known 
and given by 

I \ sinhz ,^ ^, 

^) = / , = • (4-4) 

V sinh^ z + e-4^ 

Consistency requires that 

m^ = m{K,R^). (4.5) 

Upon combining (14. 3 P with (14. 5 P one obtains an equation for Hq [or equiv- 
alently mo]. The solution is a function of the two system parameters K 
and ?7 and given by 

[ / tanh^ir-e-^^ y 
tanhi7o(^,^) = < V l-e-4^ / ' " (4.6) 

[ 0, ir<i^c, 

in which there appears a critical coupling = J /T^ that is the solution o^ 

tanhr^Kc = e-2^^ (4.7) 

The magnetization mQ{K,r)) follows directly from (14. 3 p and (14. 6p . For T ^ 
TjT it vanishes as ttiq oc (Tc — TY with a classical exponent /3 = |. For later 
use it is worthwhile to notice that also Hq{T) oc (T — Tc)^/^ when T < Tc . 

The DB property found below equation ( 14. ip now allows us to conclude 
that for — )■ cxD the stationary state distribution Pst,o('") s) is the Boltzmann 
distribution corresponding to (14. 2p . that is, 

Pst,o(r,s)=Aroe-^"('^'^)/^ (4.8) 

where Ao is the normalization. In reference [1] several system properties were 
calculated in this N ^ oo limit by averaging with respect to Pst,o('^; -s). 



5 Expansion procedure for P^t^r^s) 

As has become clear in section [3l the inverse system size 1/A is a measure of 
the degree of DB violation. In the present case this will lead us to attempt 
to find the finite N stationary state by expanding around the known N = oo 
solution (14. 8p . which will play the role of the zeroth order result. At the basis 

^Equation (|4.7p may be rewritten as sinh(2 Ji/Tc) sinh(2 J2/rc) = 1, which shows, as 
was also noticed in reference [1], that Tc is exactly (but accidentally) equal to the critical 
temperature of Onsager's square Ising model with horizontal and vertical couplings Ji and 
J2. 
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of the expansion is the hypothesis, to be confirmed self-consistently, that the 
fluctuations of the chain magnetizations, defined by 

Sfi{r) = jj,{r) — niQ , 6fi{s) = — mp . (5.1) 

are of order N~^^'^. 

A naive attempt to set up the expansion would be to notice that the 
transition probability fl2.7ap can be written as a sum of its average and a 
correction, Wj{r; s) = Wjfi{r) + Wj{r; s), where Wj^ir) is given by (14. ip and 
Wj{r]s) = wf{r) X {—^rj)Sjj{s)tanhr]K is of order N~^^'^. One might then 
think that there exists a corresponding expansion Pst(r, s) = Pstfiif^, s) [1+. . .]. 
However, the dot terms turn out to be of order 0{1) as — ?• oo, which is a 
sign that this is not the right way to expand. The reason for this failure is 
that Pst is the exponential of the extensive quantity "Ho; one should therefore 
ask first if this exponential contains any corrections of less divergent order in 

before attempting to multiply it by a series of type [1 + . . .]. In the next 
section we describe how the expansion can be set up successfully. 

Knowing how to calculate higher order corrections to the stationary state 
distribution, although certainly of diminishing practical interest, has a def- 
inite theoretical merit. What we will find in the end is that in fact to first 
order in the expansion detailed balancing continues to hold, but with respect 
to a Hamiltonian T-L^^^{r,s) that acquires a first order correction. In section 
Owe present the solution, to be denoted as Pst\f^, s), of the stationary state 
to first order. In section [7] we will show how higher orders can be calculated 
and find that from the second order on DB violation appears. Section [7] also 
provides the demonstration of the correctness of the expansion. 



6 Stationary state Pst(r, s) to first order 

We use the upper index '(1)' to indicate any quantity correct up to first order 
in the expansion. We will prove that the correct expansion takes the form 

Pst(r, s) = P,P(r, s) [1 + gi(r, s) + q^ir, s) + . . .] , (6.1) 

where the qt (/c = 1, 2, . . .), that we will show how to determine later, are of of 
order 0{N~^f'^) and where P^l\r,s), which includes a first order correction 
to the zeroth order result, is explicitly given by 

i^?(r,.)=Ar«exp(-^^ii^), (6.2a) 
f = — f goN6fi{r)6fi{s), (6.2b) 
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(^0 = cosh^ i^o tanh?7-ft', (6.2c) 

in which A/"*-^-* is the appropriate normahzation. The second term on the RHS 
of fl6.2bp is a correction to the zeroth order effective Hamiltonian. It is 0{1) 
ioY N ^ oo and, since it is proportional to go, it vanishes as expected when 

7] = 0. 

In order to demonstrate fl6.ip -f ET^ we split W according to 

oo 

>V = >V« + ^Wfc, (6.3) 

k=2 

where we take for W^^-* the matrix with the factorizing transition probabilities 
that ensure detailed balancing with respect to 'H^^\ and in which the Wfc will 
be defined shortly. Expression (16. 2p for Hamiltonian H*-^^ shows that a spin 
rj is subject to a total field Hq + go6fi{s). Hence by analogy to (14. ip the 
transition probabilities that enter W*-^^ are 

wf\r; s) = wf{r) x i[l - rj tanh{i7o + go^f^is)}]. (6.4) 

We then have by construction that 

>V(i)P3^) = 0, (6.5) 

which is the combined zeroth and first order result. It may be obtained in 
explicit form from (13. 2p by the substitutions Wj t-)- wj^^ and Pst ^ P^t^ ■ 

A remark on terminology is in place at this point. Since the zeroth and 
first order will often be combined, we will refer to equation (16. 5p as describing 
the 'leading order'. The terms qi, q2, ■ ■ ■ in the series (16. ip will be referred to 
as 'higher order' corrections. 



7 Stationary state to higher orders 

The validity of the expansion procedure of this section hinges on our being 
able to show that the corrections take effectively the form of the series of 
in (16. ip . where the terms are proportional to increasing powers of iV~^/^. 

7.1 The perturbation series for Pst{r,s) 

In order to show that the higher order corrections to Pst can be expressed as 
the series of equation (16. ip . we must first define the Wk in equation (16. 3p . 
Let us define 6wj{r; s) by 

Wj{r; s) = wf'\r; s) + 6wj{r; s). (7.1) 
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Starting from (17. ip we employ the explicit expressions f l2.7ap and (16.40 for 
Wj and wfK respectively, perform a straightforward Taylor expansion in 
and still use f l4.3p to eliminate rriQ in favor of Hq. This leads to 

5wj{r; s) = wf{r) x |[1 — rjfi{s) tanhr/K] — |[1 — tanhjifo + go^K-^)}] 
= w'^{r) X {—^Tj) ^{niQ + tanhr^i^T — tanh{ifo + doSfiis)} 

oo 

= wf{r) X (-ir,)^afc5/(s) 

k=2 

oo 
k=2 

where the last equality, supposed to hold term by term in k, defines Wj^k and 
shows that it is of order A^~'^/^. In the third line of (17. 2p the vanishing of the 
term linear in 6fi has of course been pre-arranged. The first two coefficients 
ttk in that line are given by 

0-2 = doi^ — tanh^ Hq) tanh Hq , 

ag = 1^3(1 -tanh2i7o)(l -Stanh^i/o)- (7.3) 

It becomes clear now that there is a qualitative difference between the high 
temperature regime T > Tc where we have Hq = 0,a2 = 0, and 

a3 = ^taiih^r]K, T>T,, (7.4) 

and the low temperature regime T < where Hq > 0, 02 > 0. 

We define the matrices in expansion (16.30 in terms of the Wj^k by 
analogy to (12. lip . Hence for T > Tc we have that W2 = 0. 



7.2 The higher order equations 

The leading order equation (16. 5 p being satisfied, we now turn to the higher 
orders. Substitution of (16. 3 p in (13. ip and use of (16. 5p leads to an expansion 
of which the first term is 

W«Pi'V + W2n? = 0, T<Te. (7.5a) 

In the high temperature phase the fact that VV2 = implies that qi=0 and 
therefore (I7.5ap is replaced by the next term in the expansion, 

^^'^Plt\2 + mP^P =0, T > T, . (7.5b) 

Either will be referred to as the 'next-to- leading order' equation. One ob- 
tains all higher-order equations in explicit form by inserting in the full sta- 
tionary state equation (13. 2p the expansions (16. ip for Pst{r,s) and (I7.ip -( rr2|) 
for Wj{r; s). 
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7.3 Equation for T < 

By the procedure indicated above we obtain for the next-to-leading order 
equation fl7.5ap the exphcit form 

j 

+ Wj^2{s^; r)P^^\r, s^) - w^- 2(3; r)P^P{r, s) 

+ wf\r^; s)P^P{r^, s)qi{r^ , s) - wf\r; s)P^l\r, s)qi{r, s) 

+ wf\s^; r)P^P{r, s^)qi{r, s^) - wf\s; r)P^^\r, s)qi{r, s) . 

(7.6) 

We wish to divide (17. 6p by P^l\r, s) and therefore have to compute 



We easily find 

2Rj{r;s) = [H^^\r\ s) - n^^\r, s)]/T 

= [Hoir^ s) - Ho(r, s)]/T - goN6fi{s)Mr^) - 5/i(r)] 
= -2Kr,{rj_, + rj^,) -2r,{Ho + go5fiis)}, (7.8) 
where we used (16. 2 p and (14. 2p . Detailed balancing says that 

w\'\r^; s)P^l\r\ s)(r^ s) = wf\r- s)P^l\r, s). (7.9) 

Using (17. Sp in the first two lines and ( I7.9P in the last two lines of (17. 6p we 
obtain 

= 5^ [^^■,2(r^;s)e-2«^('-''^)-«;,,2(r;s) 
j 

+ Wj,2{s^] r)e"2^^ - Wj^2{s] r) 
-f- wf\r; s){qi{r\s) - qi{r, s)} 

+ wf\r;s){qi{r,s^)-qi{r,s)} . (7.10) 

The expression in the first line of (I7.10p may be rewritten as 
Wj^2{r^; s)e~'^^'^'"'''^ - Wj,2(r; s) 

(7.11) 

= wf{r) X i[l — rj tanhifo] x 4:rj6fi^{s)gQ tanhi^o , 
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of which the first two factors on the RHS are again exactly wj^o . In f lT.lOp 
Wj^^ is of order but contains corrections of higher order in N~^^'^. In 
f l7.10p . to leading order in A^~^/^, we may therefore replace it by its — cxo 
limit, that is, by wj^ defined by (14. ip . When we substitute (17. lip in f l7.10p 
and apply to Wj^^ the — ?■ oo limit, we obtain the final form of the equations 
for the next-to-leading order correction to the stationary state, 

= |^Wj^o('^) X 4rj5fQ tanhiJo (^/^^(s) 

j 

+ Wjfl{s) X 4sj(?o tanh Hq Sfj,'^{r) 

+ U!j,oi^){'liir\ s) - qiir, s)} 

+ ^/;,-o(s){gi(r,s^)-gi(r,s)}]. (7.12) 

7.4 Equation for T > Tc 

For T > Tc we have 02 = whence gi = 0. Equation (17.5bp . when rendered 
explicit, leads to expressions that are identical to successively (17. 6p . (I7.10p . 
and (I7.12P apart from the substitutions gi 1— )■ q2 and Wj^2 ^ wj^^. In this 
case Wjfl{r) = |w^(r) where wf{r) is given by (12.50 . and 6fi = fi since 
Hq = mo = 0. Hence instead of (I7.12p we get 



+ wf{s) X 4sja3fi^{r) 

+ wf{r){q2ir\ s) - Qiir, s)} 

+ wf{s){q2{r,s^)-q2{r,s)}'^. (7.13) 

Finding the solutions of (I7.12p and (I7.13P will be the subject of the next two 
subsections. We will first consider the easier case of T > and then the 
case T < Tc. 

7.5 Solution for T > 

We start with the high temperature phase, where equation (17.130 applies. 
Detailed balancing would be satisfied if the expression under the sum on j 
were zero, that is, if we had 

q2ir^ , s) - q2ir, s) = -a^rj^i^i-s), 

q2{r,s^) - q2{r,s) = -asSjfi^^r). (7.14) 
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It can easily be shown that it is impossible to satisfy these equations. How- 
ever, they suggest that we look for a solution q2 of the form 

g2(r, s) = NC2a3[fi{r)n%s) + n{s)fi%r)] (7.15) 

where only the constant C2 is still adjustable. The difference q2{r^ , s)—q2{r, s) 
is easy to calculate, but we are interested only in its leading order. This leads 

to 

q2ir^,s)-q2ir,s) = -2C2asr,[fi^is) + 3fiis)fi\r)] + OiN-^), (7.16a) 
q2{r,s')-q2{r,s) = -2C2asS,[fi%r) + 3fi{r)fi\s)] + 0{N-^). (7.16b) 

It should be noted that whereas fl7.15p is of order A^~^, the differences (17.161) 
are of order N~^/'^ . We now need 

J^w^f (r){g2(r^s)-g2(r,s)} = -2C2a^ {^^firy^j [^?{s)+?>^i{s)^i\r)]. 

(7.17) 

With the aid of the explicit expression for wfij-) one evaluates easily 

Y,wf{r)r, = \{l-^)N^^{r). (7.18) 
j 

we see that the equation is satisfied for €2 = \- Hence from (17.15P we get 

92(r,s) = ^^N{tiin\i^ T]K)[ii{r)ii^{s) + ii{s)ir^{r)]. (7.19) 
This is of order A^~^. 



7.6 Solution for T <% 

In the low-temperature regime equation fl7.12p applies. In order to solve this 
equation we now postulate 

gi(r, s) = NCM^f^Hr) + Sfx^s)] (7.20) 

where Ci is an adjustable constant and 

62 = 4a2/(l - tanh^ Hq) = Agl tanh Hq . (7.21) 

Expression f l7.20p is of order A^~^/^. Instead of (17.160 we now have the dif- 
ference 

gi(r^ s) - gi(r, s) = -6Cib2rj6iJ,\r) + 0{N-^/'^). (7.22) 
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which is of order ^. The first two hnes of f l7.12p require that we evaluate 
J2 ^j'Oi'^yj = lJ2^^^ - '^i+i)] [1 - tanh Ho]rj (7.23) 

3 3 

Unhke the sum in f l7.17p . this is not a sum of zero-average random terms. It 
will produce a result of order A^, which we may replace by its average. This 
yields 

"^Wjfiiryj = iA^[(l -7)mo - (1 -7aH)tanhifo] 

3 

= NG, (7.24) 

where the last equality defines G and where an is the nearest neighbor spin- 
spin correlation (rjrj^i) of a ID Ising chain in a field as described by "Hq 
[equation fl4.2p ]. Expression fl7.24p . contrary to its T > Tc counterpart fl7.18p . 
has no spin dependence and is therefore equal for the r- and s- spins. The 
first two lines of fl7.12p . to be denoted Si, become 

Si = ANGg^ tanh Hq [djj^s) + 6iJ,\r)] . (7.25) 

We use f l7.22p to write the last two lines of f l7.12p as 



5*2 = — 6C162 



= -QNGCib2[6i2\r) + 6i2\s)]. (7.26) 

The stationary state equation f l7.12p may the be written as 5*1 + 5*2 = and 
we see that it is satisfied for Ci = |. 

gi(r,s) = fAr^^tanhi^o [6 jj^r) + 6 fi^s)] . (7.27) 
7.7 Section summary 

We have studied in the preceding subsections the large- expansion of the 
stationary state distribution Pst('", s) of the infinite velocity CRIC defined in 
section O We have shown, for T < Tc and T >Tc separately, the existence of 
a series of correction terms qk that multiplies the leading order result P^l"^ in 
(16. ip . which itself is again composed of a zeroth and a first order contribution. 
This expansion also furnishes the necessary proof that the prefactor P^l^ 
represents indeed the 'leading order' behavior. We have determined explicitly 
the first nonzero correction term in this series: qi for T > Tc and q2 for 
T < Tc 
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When looking ahead beyond this leading order correction, it appears that 
the Qk (for k >2 when T < Tc and for k > 3 when T >Tc) involve not only 
5/i(r) and 5/i(s), but also energy fluctuations such as A^^"*^ ~ ^h), 

if not longer-range correlations. Therefore, even though on the basis of the 
results of this section one might be tempted to postulate a general solution 
of the simple type Pst('^, s) = Pst\f, s)Q{6iJ,{r), 6jj{s)), it is unlikely that the 
true Pst('^, s) is of this form. 



8 Stationary state averages 

Stationary state averages of observables A{r, s) are averages with respect to 
Pst{r, s), so that using (16. ip and (I6.2ap we have 

= + + (8.1) 

where (. . .)*-"^^ indicates an average with weig ht P^t\r, s) [equation (^], the 
second line results from a straightforward expansion, and 

T > Tc, 

(8.2) 

T < T, , 

for the lowest order nonzero terms in the expansion. Although the qk are 
accompanied by increasing powers of N~^^'^, the order in N~^/'^ of each of the 
terms in the series (18. ip must be analyzed for each observable A separately. 



8.1 Integral representation of the partition function 

The denominator in the first line of (18. ip is a normalization factor to which 
we may refer (although slightly improperly) as the partition function Z. In 
order to find expressions for the averages (. . .)^^) in the second line of (18. ip . 
we begin by evaluating Z to leading order, 

Z^^\K,H,,g,)^Y.^-^'^'^'''^l^^ (8.3) 

r,s 

with 'H(i) given by (]6.2bp in which one should substitute (14.21) and (15. ip . To 
this order (18. 3p is a true partition function, viz. the trace of a Boltzmann 
factor. The notation Z^^\K, Hq, go) is meant to indicate that we wish to 
consider this quantity as a function of three independent parameters, ignoring 
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for the moment expression (I6.2cp for go . The r- and s-spins in (18 .Sp may be 
decoupled by the integral representation 



oo fOO 

dx 



oo J — oo 



dy e 



-g-^N(x^+y^) 



^-(x+iy)Nmo J2j rjrj+i + {Ho+x+iy) Y.j 



{x-iy)Nmo ^ ■ SjSj+i + (HQ+x-iy) s 



3 '3 



in which mo = m{K,Ho) follows from (I4.3p and (14. 6p . The two factors in 
brackets in (18. 4 p are seen to be the partition functions ({K, Hq + x ± iy) of 
independent standard Ising chains in magnetic fields Hq + x ± iy. Hence 

AT ;"00 j'OO 

Z« = — / dx / dye-^o-^^(^^+i/^)-2^^"^«|C(i^,i^o + a; + iy)|'. (8.5) 

T^90 J-oo J-oo 



We recall that 
where 



X±{K,B) = e 



K 



cosh B ± Vsinh^ B + e-^^ 



(8.6) 
(8.7) 



are the transfer matrix eigenvalues. 



8.2 Stationary point and fluctuations 

The X and y integrals in (18. 5 p are easily evaluated by the saddle point mey- 
hod, In the limit of large N, we may neglect in (18. 6p the exponentially small 
corrections due to A_ and get from (18. 5p 

ZW ^ / dx dy e-^^^-'V)^ (8.8) 

TT^'O J-oo i-oo 

where 

^(x,y) = g^\x^ + y^) + 2xmo - log | A+(ir, ifo + x + iy)\^ (8.9) 

Let {x*,y*) denote the stationary point of the integration in (18. 8p . The 
stationary point equations = ^y = ^ can be expressed as 

g^\x* ± ly*) = m{K, Ho + x* T iy*) - , (8.10) 
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with the magnetization m(K,B) = X^^{K, B)d log X+{K, B)/dB given by 
f l4.4p . For reasons of symmetry the stationary point must have y* = 0. This 
reduces flS.lOp to the single real equation 

g^'x* = m{K, Ho + x*) - m{K, Ho), (8.11) 

where we used that mo = m{K,Ho) [equation fl4.5p ]. Equation (18. lip has 
for all Ho the obvious solution a;* = 0. We investigate the stability of the 
stationary point {x*,y*) by calculating the matrix of second derivatives, 

J-:, = 2[g,' - xiK, Ho)], J-;, = 2[g,' + xiK, Ho)], 

J^xy = J^yx = 0) (8.12) 

where the asterisk indicates evaluation in the stationary point and where 
x{K, B) = dm{K, B) / OB is the magnetic susceptibility. We obtain the eigen- 
values J-"*^. and J^yy explicitly by substituting in (I8.12p for go the expressions 
(I6.2cp and for x the expression 

, cosh B 
X{K,B) = - — ^, 8.13 

(sinh' 5 + e-4^) ^ 



T>Te, 

(8.14) 



where (14. 4p has been used. This yields 
' 2{e-^^ Ttanhr]K) 



= I 

xx,yy ^ 



e 2^ tanh r]K ' 
2(1 -tanh^ r/K) (tanh^ r/iT ^ e"^^) 
(1 - e-4^) tanh^ r]K 



in which the upper (lower) sign refers to the xx (to the yy) derivative. It 
can be seen that J^yy is positive for all temperatures, but that J-"*^ , which is 
positive in both the high and the low-temperature phase, vanishes as T — )■ Tc. 
Hence for all T ^ Tc the stability is ensured by the quadratic terms in the 
expansion of T{x,y) around the stationary point. 



8.3 Free energy 

We are now in a position to calculate various physical quantities of interest. 
The first one will be the interaction free energy per spin between the two 
chains which (divided by T) will be called F^■^^t. It will turn out to have an 
expansion 

Fint = + iV- + . . . (8.15) 

To show this we pursue the calculation of Z^^^ begun in (18. 8p . We there 
substitute the expansion 

Hx, y)=r + \7l,x^ + i + . . . . (8.16) 
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We can then carry out the integrations in (18. 8 p by the saddle point method 
and find that only the quadratic terms in (18.161) contribute. The result has 
the form 

^(1) _ g-A^^*-/int [1 + 0{N-^)] (8.17) 



where 
and 



ATjr* = NT{0,0) = -2N\og\+{K,Ho) {8.U 



fUK,v) = I log [l-glx\K,Ho)] . (8.19) 

Here J-"* is the free energy (divided by T) of two independent Ising chains in 
an effective field Hq. Since Hq is proportional to the coupling rjK between the 
chains, the field dependent part of J-"* actually represents the bulk interaction 

int 



free energy NF^^^ between the chains, that is 



0, T > Tc , 



and furthermore fi^t{K, 7]) is a residual interaction free energy between them 
which remains of order iV° as — > oo. The energy that one drives from it 
has a cusp singularity and hence the exponent a = [T]. 

Beyond this leading order result we obtain /int explicitly in terms of the 
two system parameters K and 7] by substituting in (I8.19P the expressions 
for go and x given in (I6.2cp and (I8.13p . respectively, and (when T < Tc) 
eliminating Hq. The result is that 

J I log (1 -e^^tanh^ r/fsT) , T > , 
^'"^"^'^^ = \ llog(l-e-«^tanh-%ir), T < . ^''''^ 

In view of (I8.20p we see that Fjnt has a linear cusp at T = Tc, and (I8.2ip 
shows that /int diverges logarithmically for T — )• Tc. In spite of this weak 
divergence, the finite size correction /int to the interaction free energy Tint 
also conforms the classical specific heat exponent a = 0. 



8.4 Finite size scaling of the free energy near Tc 

We will show how our approach allows for finding the finite size scaling 
functions. By the way of an example we consider the singular part of the 
free energy. For T ^ the quantity /int diverges due to the second order 
derivative J-"*^ becoming zero. In order for the integral (18. 8 p combined with 
(I8.16P to converge at T = Tc, we have to include higher order terms in the 
expansion (I8.16p . We will write 

•^i^yU) — J' ~\~ 2*^xx'^ ~^ ~^ '%^xxx^ ^ '2i^ xxxx^ + . . . (8.22) 
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and will argue below that near Tc the terms not exhibited explicitly in this 
series are of higher ordeij^. In order to find the coefficients in fl8.22l we perform 
a straightforward derivation of (18. 9p and set x* = y* = Q. We then define 

e = = ^ (8.23) 

which, in the vicinity of , leads to 

Ho = B,e'/^ + 0{e) (8.24) 

where from (14. 6 p we have 

B = < (8 25) 

" ~ \ 2e-^^-{ri + l/smh2K^)K^, T < . 

When using (I8.24p in the coefficients found above we obtain 

J-:, = a±e + 0(e2), 

K... = c + Oie), (8.26) 



where 



a± = <! ^ |(e^^'^-l)(r/ + l/sinh2ir,)ire, ^ ^ ' 



0, T > Tc , 

4(3e^^^ + 1/ sinh2is:c)i^c , T <T,. 

c = 6e2^% (8.27) 

We substitute the explicit expressions (I8.26P in (I8.22p and use that expansion 
in the integral (18. 8p . When we introduce the scaled variables of integration 
u and V defined by 

X = N-^/\ y = N-^/^v, (8.28) 
as well as the scaling variable 

r = eA^^/^ (8.29) 



"Terms with an odd number of y derivations vanish by symmetry. 
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the factor N disappears from the exponential. After carrying out the Gaus- 
sian integration on v we get 

^(i)^e-^^*x^^— ^Z(r), (8.30) 



. 27r 

vahd in the scahng hmit TV — ?• oo, T — )■ with r fixed, and where Z is the 
scahng function 

/oo 
du exp [~la±\T\u^ - lb±{-rY/^u^ - j^cu^] . (8.31) 
-oo 

It is of a type that occurs standardly in problems with mean field type critical 
behavior; they have been studied recently by Griineberg and Hucht [8]. It 
has the limiting behavior 

Z{0) = J^^due-^^'/^\ r^O, 

Z(r)^<; y/2 (8.32) 

Vl I , r — )■ ±oo. 



Upon combining f l8.17p and f l8.30p we find that 

/i,t(i^, r^) = -i log iV - log Z{tN'/^) - i log + . . . , (8.33) 

again valid in the scaling limit, and where the dots stand for terms that 
vanish as — t- oo. It follows, in particular, that equation f l8.2ip may now 
be completed by 

/i,t(i^c,r7)~ilogAr + logZ(0) + ..., T = Te, AT ^ oo, (8.34) 
where the dots stand for terms that vanish as A^ — ?■ oo. 



8.5 Susceptibilities 

Of primary interest are the correlations between the fluctuations of the mag- 
netizations in the two chains. We set as before 6fi = fi — uiq. The general 
expression that we will study here is 

= {6fi''{r)6i2'{s))^^^ + ... (8.35) 

where the dots in the last line, obtained according to (18. ip . represent higher 
order terms. Special cases that we will consider are the cross-chain suscepti- 
bility Xint and the single-chain susceptibility Xsin, defined as 

Xint = Nxii = N{5fxir)5fx)is)), (8.36a) 
Xsin = Nx2o = N{6fi\r)), (8.36b) 
in which, of course, the latter is also equal to Xo2 by symmetry. 
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8.5.1 Cross-susceptibility 



We first consider the correlations between the fluctuating magnetizations of 
the two chains. The cross-susceptibility Xint is the quantity most character- 
istic of these correlations. From equations (I6.2bp and f l8.3p it is clear that 
Xint = c^log Z^^ydgo where the derivative has to be evaluated at flxed K and 
Hq, considering qq as an independent parameter in (18.41) . Doing the calcula- 
tion for Z(^) given by flHTTD . f l838|) . and f lHTTOj) . we observe that = J'(0, 0) 
is independent of go so that 



XintiK, ri) 



dfint 

dgo 



9oX 



9oX^ 



tanh rjK 



e-^^ - tanh^ r]K 
e-8^fl - 



tanh^ r]K) 



(tanh e'^^) (tanh^ rjK - e-^^] 



T>T,, 
(8.37) 

T<Tr. 



For T — )• Tc this quantity diverges as |T — Tc|~^'"* with 7int = 1. It is a signal 
that at T = Tc this correlation scales with another power of A^. A scaling 
function for Xint may be derived from the one for /int , but we will not try to 
be exhaustive. 

Since at speed f = cxd all index pairs (i, j) are equivalent, the correlations 
between the r- and the s-spins are given by 



(8.38) 



8.5.2 Single-chain chain susceptibility 

The single-chain susceptibilities Xsin is deflned in equation (18.361) . Let us 
now consider the general expression (18.351) for Xki , for which the appropriate 
approach differs slightly from that of the preceding subsection. One may 
generate insertions Sfj,''{r [or 5/i^(s)] in the integral (18.51) by passing from 
X and y to the two independent variables z = {x + iy) and z = {x — iy) 
and letting N-^d^/dz^ [or N-^d^/dz^] act on q-^^^'^'' Z {K , Hq + z) [or on 
e~'^^^'^° Z {K , Hq + z)]. We flnd, using (18. 6p and neglecting again the effect 
of A_ which is exponentially small in A^, 

^''tTT [e-^^'"»Z(ir, Ho + z)] = Uz) Z{K, Ho + z), (8.39) 
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in which 



^0 


[z) 


= 1, 




Jl 


[z) 


= m — mo , 




■h 


[z) 


= {rh- mof + N~^x-, 




J3 




= (m - mo)^ + 3iV~^(m 


- mo)x + N-^x', 


J4 




= (m — mo)^ + 6A^~^(m 








+3N~'x' + N-^x", 





mo)x' 



]A0) 



where, in this formula, we abbreviated m = m{K, Hq + z) and x = x{K, Hq + 
z) [see equations ( I4.4p and ( 18.13P ] in order to emphasize the z dependence 
of these quantities, and where the primes on x stand for differentiations 
with respect to Hq. Equations f l8.39p and (18.401) of course have counterparts 
obtained by letting r h-)- s, £ and z ^ z. When (I8.39P is substituted in 
(I8.35P we obtain 

XM = {Jk{z)JfX-z))'^'^ + ..., (8.41) 

where the dots stand for higher-than-leading order terms in the A^~^ expan- 
sion. 



By virtue of equations (I8.4ip and (I8.40p it follows that 
X20 



= {{m{K,Ho + z)-moY) + N-\x{K,Ho + z)) (8.42) 

We now expand m and x for small z anticipating that upon integration with 
weight exp(— A^J-") each factor z"^ will, to leading order, produce a factor A^~^. 
After multiplication by this yields 



Xsin(i^, V) = Nx'{K, Ho) + x{K, Ho) + 0{N-' 



(8.43) 



Anticipating again that each factor z or z will produce a factor A^~^/^, we 
see that all terms exhibited explicitly on the right hand sides in (18.480 are 
of order A^~^. We have replaced the averages (. . .)^^\ which are with respect 
to exp{—NJ'{x,y), by averages (. . in which F{x,y) of equation (18.90 is 
replaced with the Gaussian terms in its expansion, shown in (18.160 . 
Upon using in (18.420 the explicit evaluations 

















six 








iv(i - ^oV) ' 



.44) 
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we arrive at 

XsUK,r])= ^ , T^Te, (8.45) 

valid in the limit N ^ oo. Hence the in-chain susceptibility Xsin is equal to 
the susceptibility of the ID Ising model enhanced by a factor (1 — g^x"^)'^ 
due to the presence of the other chain. 

Using expressions f l6.2cp and fl8.13p for and X; respectively, we may 
render fl8.45p explicit in terms of K and rj and get 



Xsm{K,ri) 



T>T, 



C ) 

r.... ._._.2„..^ (8-46) 



e-^^ - tanh^ r]K 
e-^^{tanhr]K){l - tanh^ r]K) 
(1 - e-4^)(tanh^ r/i^ - e'S^) 



For T — )• Tc the susceptibility Xsin diverge as (T — T^)^^' with, again, the 
classical critical exponent 7 = 1. For 77 = (whence = 0) the first one of 
equations (18.460 reduces to the standard susceptibility of the zero field ID 
Ising chain. 

In agreement with the symmetry of the problem, Xint is odd and Xsin is 
even in rj. Both above and below Tc one easily verifies that in agreement 
with Schwarz's inequality we have Xint/Xsin < 1- 



8.6 Spontaneous magnetization 

For T > Tc symmetry dictates that the magnetization (yu(r)) and (/i(s)) 
are zero to all orders. However, for T < Tc the magnetization fi{r) = 
^jLi '^j has, to leading order, a Gaussian probability distribution of 
width N'^/"^ around mo(K, Hq). As a consequence (5/i(r)) vanishes to order 
A^~^/^. However, to order there appear nonzero corrections terms to 
(/x(r)). As an application of equation (18. ip we calculate in this subsection 
these correction terms. 

Upon using (18. ip for the spacial case A = 6fi{r) and inserting in it the 
explicit expression (I7.27P for qi we obtain 

(5/i(r)) = ((5/^(r))« + iNg'.t^nhHo [{Mz)f^ + {Ji{z)Mz)f^] . (8.47) 

When substituting (I8.40p in the second term of (I8.47P we see that we need 

(Mz))^'^ = x^^')g^ + 6Ar-ix^(z2)«+3Ar-V + 0(iV"^/^), 
(Ji(^)J3(z))« = x'{zz')^i^ + SN~\^{zz)'^^^ + 0{N''/'). (8.48) 

We have replaced the averages (. . .)^^^ by averages (. . for the same reasons 
as in the preceding subsection. Taking into account again that each factor 
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zoT z brings in a power iV~^/^, we see that all terms explicitly exhibited on 
the right hand sides of equations (18.481) are of the same order in A^, namely 
(9(A^~^). The Gaussian averages are easily calculated and we are led to 

(Mz))^'^ + {JiizUsiz)^'^ = + 0{N-'/'). (8.49) 

^ (1 - 9oX ) 

We should now evaluate the first term on the right hand side of (I8.47P , namely 

{6f^{r))^'^ = {Mz))^'^ = x{z)^''^. (8.50) 

The Gaussian average {z)q^ vanishes on account of symmetry. However, 
when the third order terms in the Taylor expansion (I8.16P of J^{x, y) are 
kept and we expand these we get after a straightforward calculation that we 
will not reproduce here, 

(5/.(r))« = |iVxx'[(a:^)S^-3(xV)S^] 

The final result for {6fj,{r)) is obtained by substitution of (I8.5ip and (18.491) 
in (I8.47p . We see that ((5/^(r)) has two contributions of order N~^. The 
contribution (Sjj^r))^^^ comes from the effective leading order Hamiltonian 
H^^^ The second contribution accompanies the violation of detailed balanc- 
ing symmetry and is therefore essentially a non-thermodynamic effect. 

8.7 Pair correlation function 

It is of interest to study the pair correlation 

g^{e) = {r,r,+,) (8.52) 

in a single chain. To that end we consider again expansion (18.11) . now with 
A = rjTj^i. Its first term may be written 

^?«(£) = Z«/Z(i) (8.53) 

where Z^^'' is given by (18. 4p but with an insertion rjTj^e in the sum on r. 
Equivalently, Z^^^ is given by the same integral as (18. 8 p but with an insertion 
g^j^\£; K,Ho + z), this quantity being the pair correlation of the ID Ising 
chain in a field Hq + x + iy. Evaluation by means of the standard transfer 
matrix method yields 



~9^^\i; K, H, + z) = m\K, H, + z)+ '"t^'^f , (8-54) 

smh [Ho + z) + e~*-^ 
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well-known in the case z = 0, in which we defined A = A_/A+ , where the 
tilde serves as a reminder of the z dependence, and where contributions 
exponentially small in have again been neglected. In order to obtain the 
desired physical correlation function g^ipj of this system we now have to 
average fl8.54p with an appropriately normalized weight exp [ — NT{x^ y)\ . 

We will consider this quantity in the high-temperature regime T > 
where = 0. Knowing that z is of order iV~^/^ we expand (18.541) for small 
z, which gives 

~g'}i\i;K,Ho + z) = 

e^^z^ ^ (tanh^ir)(l - e^^z^) (_(e-4^ + e^^i)z^) + 0{(^:^) 

(8.56) 

To leading order the average on z may be carried out with the weight 
exp [ — NF{x, y)\ in which the expansion F is limited to its quadratic terms. 
Straightforward calculation yields 

2 3 

^jv(^) =tanh^K+(l-tanh^i^) N-^+OiN-^), T > T, , (8.57) 

valid for N oo at fixed £, where as before x stands for the susceptibility 
x{K, 0) = e^^ of the ID Ising chain and where go = tanhriK. In the scaling 
limit £, — 7- oo with a fixed ratio one obtains 

2 3 

gN(i) ^ ( tanh^ K)(l)(iN-^) + , A^-^ T > , AT ^ oo, (8.58) 

1 - 9oX 

in which each of the two terms is valid up to corrections of relative order 
A^~^ and in which is the scaling function defined by 

<P'i^) = -, 2~2ff 2 2 • (8-59) 

We observe the noncommutativity 

N/2 oo 

lim Yl 9N{i) ^ Yl 9N{i). (8.60) 

e=-.]\f/2+i e=-oD 

The right hand side of this inequality is equal to x{K, 0) whereas the right 
hand side is equal to x{K, 0) -|- Xint{K, rf). 

We conclude by noting that the pair correlation function may also be 
studied to higher order in A^~^ in the low-temperature regime. For T < 
the fiuctuations of the magnetic field z are asymmetric and greater care is 
required. We will not include such a calculation here. 
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9 Traffic model 



Motivated by an interest very different from that of references [H [2] we re- 
cently introduced a new traffic model describing vehicles that may overtake 
each other on a road with two opposite lanes p]. That work shows the ap- 
pearance of a phase transition when the traffic intensity, supposed equal on 
the two lanes, attains a critical value. Above the critical intensity the sym- 
metry between the two traffic lanes is broken: one lane has dense and slow, 
the other one dilute and fast traffic. The study of reference [1] invoked a 
mean-field-type assumption that couples the velocity of a vehicle in a given 
lane to the average of the vehicle velocities in the opposite lane. This as- 
sumption was justified by the argument that a vehicle in one lane encounters, 
in the course of time, all vehicles in the opposite lane. Although there is no 
one-to-one correspondence between the two models, they share essentially the 
same features, as may be seen as follows. For J2 < the two chains of the 
CRIC studied here have opposite spontaneous magnetizations; up-spins may 
then be regarded as the vehicles of the traffic problem; they will be denser 
in one chain (traffic lane) than in the other. The CRIC is more amenable 
to analysis than the traffic model. It was shown analytically [H [2] that the 
CRIC phase transition disappears when v is finite. Our simulations [TU] of 
the traffic model have shown, nevertheless, that this problem is close to the 
critical point v = 00. This explains the critical-point-like phenomena that 
we observed, namely ffuctuations that last longer than the simulation time. 

10 Conclusion 

We have considered in this paper the nonequilibrium steady state (NESS) of a 
model consisting of two counter-rotating interacting Ising chains introduced 
by Kadau et al. [2] and by Hucht [T]. The model is related to a road 
traffic model studied earlier by ourselves [9J. Its dynamics is governed by a 
master equation parametrized by two interaction constants J/T and rj. The 
model has a phase transition, known to be of mean field type, at a critical 
temperature T = Tc . 

Starting from the master equation we have shown that in the limiting 
case of a relative velocity f = 00 of the two chains, the stationary state dis- 
tribution Pgt may be studied in an expansion in powers of the inverse system 
size N~^. Knowing this distribution we have calculated, also as expansions in 
N~^, of averages of physical interest: the interaction free energy between the 
chains, the in-chain and cross-chain susceptibilities, the correlation function 
(for T > Tc), and the spontaneous magnetization (for T < Tc). We have 
shown how near criticality scaling functions may be explicitly calculated. 

Whereas to leading order the force exerted by one chain on the other is 
that of an effective magnetic field Hq, the expansion requires that we 
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take into account the fluctuations of tliis field around its average. It then 
appears that to leading order the dynamics obeys detailed balancing with 
respect to an effective Hamiltonian, as was found by Hucht [Ij, but that to 
higher order in the expansion the detailed balancing is violated. 

In this work we have addressed many different, albeit interrelated, aspects 
of the finite-size CRIC. We have not tried to be exhaustive and have not 
considered, for example, energy dissipation. Similarly, the parallel problem 
with open boundary conditions has been left aside. We hope that the results 
of this work will be helpful in guiding the study, which we believe to be 
worthwhile, of the finite-velocity {v < oo) version of the model. 
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